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Abstract — We consider minimal-rate sampling schemes for 
infinite streams of delayed and weighted versions of a known 
pulse shape. The minimal sampling rate for these parametric 
signals is referred to as the rate of innovation and is equal 
to the number of degrees of freedom per unit time. Although 
sampling of infinite pulse streams was treated in previous works, 
either the rate of innovation was not achieved, or the pulse shape 
was limited to Diracs. In this paper we propose a multichannel 
architecture for sampling pulse streams with arbitrary shape, 
operating at the rate of innovation. Our approach is based 
on modulating the input signal with a set of properly chosen 
waveforms, followed by a bank of integrators. This architecture is 
motivated by recent work on sub-Nyquist sampling of multiband 
signals. We show that the pulse stream can be recovered from 
the proposed minimal-rate samples using standard tools taken 
from spectral estimation in a stable way even at high rates 
of innovation. In addition, we address practical implementation 
issues, such as reduction of hardware complexity and immunity to 
failure in the sampling channels. The resulting scheme is flexible 
and exhibits better noise robustness than previous approaches. 



I. Introduction 

Digital processing has become ubiquitous, and is the most 
common way to process analog signals. Processing analog 
signals digitally must be preceded by a sampling stage, care- 
fully designed to retain the important features of the analog 
signal relevant for the processing task at hand. The well known 
Shannon-Nyquist theorem states that in order to perfectly 
reconstruct an analog signal from its samples, it must be 
sampled at the Nyquist rate, i.e., twice its highest frequency. 
This assumption is required when the only knowledge on the 
signal is that it is bandlimited. Other priors on signal structure 
ID J Ell which include subspace |[3]-l5|, sparsity ||6l-|l8l, or 
smoothness priors HI, ||9l, ifTOil . can lead to more efficient 
sampling. 

An interesting class of structured signals was suggested 
by Vetterli et al. ifTTI . ifTZl . who considered signals with a 
finite number of degrees of freedom per unit time, termed 
by the authors as signals with finite rate of innovation (FRI). 
For such models, the goal is to design a sampling scheme 
operating at the innovation rate, which is the minimal possible 
rate from which perfect recovery is possible. A special case 
that was treated in detail are signals consisting of streams 
of short pulses. Pulse streams are prevalent in applications 
such as bio-imaging 1 13|, neuronal activity and ultra-wideband 
communications. Since the pulses are highly compact in time, 
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Standard sampling methods require very high sampling rates. 
The main idea is to exploit the fact that the pulse shape is 
known, in order to characterize such signals by the time- 
delays and amplitudes of the various pulses. Targeting these 
parameters allows to reduce the sampling rate way beyond 
that dictated by the Shannon-Nyquist theorem. In fact, it was 
shown in lfT4ll that the sampling rate can also be reduced for 
short pulses with unknown shape. 

Following this parametric point of view, a sampling scheme 
for periodic streams of pulses was developed in ifTTI . lfT2ll . 
which operates at the innovation rate. It relies on the observa- 
tion that the time delays and amplitudes can be recovered from 
a set of the signal's Fourier series coefficients. This follows 
from the fact that in the frequency domain, the problem 
translates into estimating the frequencies and amplitudes of 
a sum of complex sinusoids (cisoids), a problem which has 
been treated extensively in the context of spectral estimation 

ma. 

In practical applications finite and infinite streams are usu- 
ally encountered, rather than periodic streams. For the finite 
case Gaussian IfTTI . and polynomials or exponentials reproduc- 
ing sampling kernels |FT6l|, were introduced. The approaches 
based on the first two kernels, are unstable for high rates of 
innovations [il3j. An alternative sampling scheme, based on 
a new family of time-limited filters, was presented in ifTsll . 
This approach exhibits better noise robustness than previous 
methods, and is stable even for high model orders. Exploiting 
the compact support of the sampling kernels in lfT3l . lfT6l . 
both methods were extended to the infinite case. Unfortunately, 
neither techniques achieve the minimal sampling rate, which 
is the rate of innovation, for infinite pulse streams. 

All previous methods were composed of a single sampling 
channel. Multichannel sampling schemes offer additional de- 
grees of freedom which can be utilized to achieve the rate 
of innovation for the infinite setting. In |17| a multichannel 
extension of the method in 1161 was presented. This scheme 
allows reduced sampling rate in each channel, but overall 
sampling rate similar to lfT6l . Another multichannel system, 
composed of two first-order resistor-capacitor (RC) networks, 
was proposed in fTSl]. However, this approach assumed that 
there is a single pulse per sampling period, an assumption 
which limits the method's time resolution. Two alternative 
multichannel methods, were proposed in [19| and 11201 . These 
approaches, which are based on a chain of integrators |fT9l 
and exponential filters 1201 . allow only sampling of infinite 
streams of Diracs at the rate of innovation. In addition, we 
show in simulations that these methods are unstable, especially 
for high rates of innovation. An alternative scheme, proposed 
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in II2TI . can operate at the rate of innovation for pulses with 
arbitrary shape. However, this approach constrains the delays 
to be constant in each period. To the best of our knowledge, a 
stable minimal-rate sampling scheme for infinite pulse streams, 
with arbitrary shape, is still lacking. 

Our first contribution treats finite pulse streams. We design 
a multichannel sampling system, based on oscillators, mixers 
and integrators. In each channel the signal is modulated 
by an appropriate waveform, followed by integration over a 
compact time interval. We derive conditions which guarantee 
that the output of each channel is a mixture of the Fourier 
coefficients of the signal. By properly choosing the mixing 
parameters, we show that the Fourier coefficients can be 
obtained from the samples. Once the set of Fourier coefficients 
is known, we use standard spectral estimation tools in order 
to recover the unknown times and amplitudes. As we show, 
the mixing scheme enables simple and practical generation of 
modulating waveforms. Furthermore, mixing the coefficients 
allows recovering the signal even when one or more sampling 
channels fails. 

Integration over a finite interval enables a simple extension 
to the infinite setting. Our infinite sampling approach leads 
to perfect reconstruction of the signal, while sampling at the 
rate of innovation. In addition, our scheme can accommodate 
general pulse shapes with finite length support. As we show 
in simulations, our approach exhibits better noise robustness 
compared to previous methods, and allows sampling at high 
rates of innovation. We also discuss a special case of infinite 
streams of pulses having a shift-invariant (SI) structure, a 
model presented in |21 1, and compare our method with the one 
in ETI . Finally, we describe how to practically generate the 
modulating waveforms and derive conditions on these wave- 
forms which guarantee perfect reconstruction of the signal. 

The scheme derived in this work follows ideas of a recently 
proposed sampling methodology for structured analog signals, 
termed Xampling |j22|, ll23l . |fT4l |. This framework utilizes 
the signal model in order to reduce the sampling rate below 
the Nyquist rate. A pioneer sub-Nyquist system for multiband 
signals [7|, referred to as the modulated wideband converter 
(MWC), was proposed in ll24l . Although treating a different 
signal model, our modulation scheme is based on concepts 
presented in [241. Both works share a similar analog front- 
end, so that the hardware prototype of the MWC, designed in 
II22I . can also be used to implement our method. 

The remainder of this paper is organized as follows. In 
Section |ll] we derive a multichannel scheme for finite pulse 
streams. Section |III] extends our results to the infinite case. 
We discuss the generation of the modulating waveforms in 
Section |IV] and present a practical sampling scheme which 
can be implemented in hardware. In Section |V] we discuss 
in more details the relations of our results to previous work. 
Numerical experiments are described in Section [VT] 

II. Finite Streams of Pulses 

A. Problem Formulation 

Throughout the paper, we denote matrices and vectors by 
bold font, with lowercase letters corresponding to vectors and 



uppercase letters to matrices. The nth element of a vector a 
is written as a„, and A^j denotes the ijth element of a matrix 

A. Superscripts (•)*, and (■)^ represent complex conju- 
gation, transposition and conjugate transposition, respectively. 
The Moore-Penrose pseudo-inverse of a matrix A is written 
as A^. We denote by diag(a) a diagonal matrix having the 
elements of the vector a on its diagonal. The continuous- 
time Fourier transform (CTFT) of a continuous-time signal 
X [t) e L2 is defined by X (uj) = J^^ x (t) e^^'^^dt. 

Consider the finite stream of pulses 

L 

x{t) =^aih{t~ti), tiel c[0,T),aieC, (1) 
1=1 

where h{t) is a known pulse shape, {ti,ai}fL^ are the un- 
known delays and amplitudes, and / is a continuous-time 
interval in [0,r). The pulse can be arbitrary as long as 

h{t-ti) =0,yt^[0,T) l = l...L, (2) 

i.e., the signal x{t) is confined to the time-window [0,T). 
This condition suggests that the support of the pulse has 
to be finite and smaller than T. In such cases the effective 
Nyquist rate will be quite large since x{t) will have a large 
bandwidth. However, the sampling rate can be reduced below 
the Nyquist rate, by noticing that x{t) is uniquely defined 
by the delays and amplitudes. Since x{t) has 2L degrees of 
freedom, {i;, ai}/|^j^, at least 2L samples are required in order 
to represent the signal. Our goal is to design a sampling and 
reconstruction method which perfectly recovers x{t) from this 
minimal number of samples. 

B. Relation to Model-Based Complex Sinusoids Estimation 

Our sampling problem can be related to the well known 
problem of a model-based complex sinusoids (cisoids) param- 
eter estimation. This approach was originally taken by Hou 
and Wu Ii25j, who were the first to show that time delay 
estimation can be converted into a frequency estimation of a 
sum of cisoids [15|. This follows from noticing that delays 
in the time domain are converted into modulations in the 
frequency domain. However, their method relied on Nyquist 
rate sampling of the signal, and their derivations were only 
approximate. Vetterli et al. 1 11 1 addressed this problem from 
an efficient sampling point of view, and derived a low-rate 
sampling and reconstruction scheme for periodic streams of 
Diracs. Their method was based on the same fundamental rela- 
tion between the delays in time and modulations in frequency. 
Following a similar path, we show that once a set of Fourier 
coefficients of the signal are known, the delays can be retrieved 
using sinusoidal estimation methods. We then design low-rate 
sampling schemes for obtaining the Fourier coefficients. 

Since x{t) is confined to the interval t G [0,r), it can be 
expressed by its Fourier series 

xit) =J2x[ky^''\ t e [0,T) (3) 

fcez 

where 

1 f 

X[k] = -j^ x^e-JT'^tdi. (4) 
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Substituting ([T]) into (|4]i we obtain 



X[k] 



L T 

ai / h{t~ti)e' 
Jo 



y 5Z 

1=1 

1 ^ 



aie 



aie 



h{t)t 



dt 



(5) 



1=1 



where the second equaUty stems from the condition in dU, and 
H{uj) denotes the CTFT of h{t). 

Denote by /C a set of K consecutive indices for which 
H {^k^ ^ 0, Vfc € /C. We require that such a set exists, 
which is usually the case for short time-support pulses h{t). 
Denote by H the K x K diagonal matrix with A;th entry 
C^k), and by V(t) the K x L matrix with kith element 
Q-jT-kt, ^ Yvjjgj-g t — {ti, . . . , is the vector of the unknown 
delays. In addition denote by a the length-L vector whose fth 
element is a/, and by x the length-A' vector whose kth element 
is X[k]. We may now write Q in matrix form as 



X = HV(t)a. 



(6) 



The matrix H is invertible by construction, and therefore we 
can define y — H^^x, which satisfies 



V(t)a. 



(7) 



Addressing the kth element of the vector y in (I7]i directly, we 
obtain 



Yk = ^ 

1=1 



(8) 



Given the vector x, d?) conforms with the standard problem 
of finding the frequencies and amplitudes of a sum of L 
cisoids. The time-delays can be estimated using nonlinear 
techniques, e.g., the annihilating filter [111 , matrix-pencil ||26J , 
Kumaresan and Tufts method [271 or ESPRIT lH (see CBl 
for a review of this topic), as long as if > 2L and the time- 
delays are distinct, i.e., ii ^ tj for all i ^ j. Once the 
time-delays are known, the linear set of equations (|7]l may 
be solved via least-squares for the unknown amplitudes. Due 
to the Vandermonde form of V(t), it is left invertible as long 
eis K > L, so that a = V^(t)y. 

C. Direct Multichannel Sampling 

As we have seen, given a vector of K > 2L Fourier 
series coefficients x, we may use standard tools from spectral 
analysis to determine the set {ti, ai}fL^. In practice, the signal 
is sampled in the time-domain, and therefore we do not have 
direct access to samples of x. Our goal now is to design a 
sampling scheme which will allow to obtain the vector x from 
time-domain samples. 

For simplicity, we set K to be an odd number, and choose 
the set IC = {- [ii'/2j , . . . , [7f/2j }. However, our results 
extend to any set /C of consecutive indices, as long as 
/C| > 2L. The Fourier coefficients X[k] can be obtained 
using the multichannel sampling scheme depicted in FigU] 



Each channel consists of modulating x{t) with a complex 
exponential, followed by an integrator over the window [0, T). 
The sample taken by the kth channel is exactly X[k], as in 
This direct sampling scheme is straightforward, and may be 
implemented using 3 basic building blocks: oscillators, mixers 
and integrators. However, from a practical point of view this 
approach has the disadvantage that it requires many oscillators, 
having frequencies which must be exact multiples of some 
common base frequency. 



►X +LfJ 



Fig. 1. Multichannel direct sampling of the Fourier series coefficients 

x[k], fc e /c. 



D. Mixing the Fourier Coefficients 

We now generalize our framework, towards a more practical 
sampling scheme by mixing several Fourier coefficients, rather 
than limiting ourselves to one coefficient per channel. The 
additional degrees of freedom offered by this extension will 
allow the design of waveforms that are easy to implement. 
Our approach is motivated by the hardware reported in |24|, 
where similar modulators are used to sample multiband signals 
at sub-Nyquist rates. 

In addition, in real-life scenarios one or more channels 
might fail, due to malfunction or noise corruption, and there- 
fore we lose the information stored in that channel. Unique 
recovery of the signal parameters from relies on having a 
consecutive set of Fourier coefficients [12|. Hence, when using 
the direct scheme, loss of Fourier coefficients, prevents us from 
recovering the signal. In contrast, when mixing the coefficients 
we distribute the information about each Fourier coefficient 
between several sampling channels. Consequently, when one 
or more channels fail, the required Fourier coefficients may 
still be recovered from the remaining operating channels. If 
their number is greater than 2L, the signal can be still perfectly 
recovered from the samples. We discuss this feature more 
thoroughly in Section HV] 

Consider a multichannel sampling scheme with p channels, 
as depicted in Fig. |2] In each channel, we modulate the signal 
using a weighted sum of cisoids given by 



-kt 



(9) 



where the weights Sik vary from channel to channel. The 
resulting sample of the ith channel is 



1 f 



'-kt 



dt 



s,kX[k]. (10) 
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x{t)- 



keic 



Fig. 2. Mixing the Fourier coefficients differently in each channel. 



To relate the samples and the Fourier coefficients, we define 
the p X K matrix S with Sik as its iA;th element, and by c the 
length-p sample vector with ith element c^. We may now write 
([Tol l in matrix form as 



c = S> 



(11) 



As long as S has full column rank, where p > K is a 
necessary condition, we can recover x from the samples by 
X = S^c. The direct sampling scheme presented earlier is a 
special case of this more general approach, with p = K and 
S = I. In Section |IV] we exploit the degrees of freedom this 
general scheme offers, and present sampling schemes which 
can simplify the hardware design, and are more robust to 
malfunctions in the sampling channels. 

We summarize this result in the following theorem. 

Theorem 1. Consider a finite stream of pulses given by 

L 

x{t) = aih{t -ti), tiel c[0,T),aieCJ = l...L, 
1=1 

where h{t) is a known pulse shape, and condition (|2]l is 
satisfied. Choose a set IC of consecutive indices for which 
H{2iTk/T) ^ 0, Vfc G K,. Consider the multichannel sampling 
scheme depicted in Fig. |2] for some choice of coefficients 
{sifel/cGK) « = 1, • ■ • Then, the signal x{t) can be perfectly 
reconstructed from the samples {ci\^^i with 

= ^ / E s^fce-^^'^Mt, (12) 

as long as p > \IC\ > 2L, and the coefficients matrix S in 
(II 11 1 is left invertible. 

As we discuss in Section IV-AI the method in ITSl can 
be viewed as a special case of Fig. |2] Since our work 
is a generalization of ifTJl . it benefits from the high noise 
robustness exhibited by ||T3l, in contrast to previous work |fTT|. 
|[T6|. It should be noted that Theorem [T] holds for a periodic 
pulse stream as well, since it can be similarly represented by 
a Fourier series, and all derivations remain intact. 

We now demonstrate several useful modulating waveforms. 

1 ) Cosine and Sine waveforms: First we set p = K. Then, 
we choose the first [K/2\ waveforms to be cos {^kt), the 
next [if/2j to be sin {^kt^, and the last to be the constant 
function 1. Clearly, these waveforms fit the form in (|9|l. It 



is easily verified that this choice yields an invertible matrix 
S. The practical advantage of the mixing scheme is already 
evident, since sine and cosine waves are real valued, whereas 
the direct multichannel scheme requires complex exponentials. 

2) Periodic Waveforms: Every periodic waveform can be 
expanded into a Fourier series. Transferring such a waveform 
through some shaping filter, e.g., a low-pass filter, we can 
reject most of the coefficients, leaving only a finite set intact. 
Consequently, such a scheme meets the form of (|9]l. In 
Section |IV] we elaborate on this concept, discuss design con- 
siderations, and show that properly chosen periodic waveforms 
yield a left invertible matrix S. 

One simple choice is periodic streams of rectangular pulses 
modulated by ±1 |24]. The strength of the mixing scheme 
over the direct one will be emphasized in Section IIV-BI We 
show that one periodic stream is sufficient for all channels, 
while each channel uses a delayed version of this common 
waveform. Therefore, the requirement for multiple oscillators 
and the need for accurate multiples of the basic frequency, 
are both removed. In addition, periodic streams are easily 
designed and implemented digitally, rather than somewhat 
complicated analog design of oscillators combined with analog 
circuits intended to create exact frequency multiples. Finally, 
if the period T changes, the analog circuit has to be modified 
substantially, whereas the flexibility of the digital design 
allows simple modifications. 

III. Infinite Streams of Pulses 
A. General Model 

We now consider an infinite stream of pulses defined by 

x{t) ^Yaih{t~ti), e R, a, e C. (13) 



We assume that there are no more than L pulses in any interval 
Im == [(m — 1)T, mT] , TO € Z. We further assume that 
within each interval condition ^ holds, and consequently, the 
intervals are independent of one another. The maximal number 
of degrees of freedom per unit time, also known as the rate of 
innovation mj, [12), is 2L/T. We now present a multichannel 
sampling and reconstruction scheme which operates at the 
minimal rate possible, i.e., the rate of innovation. 

Consider an extension of the sampling scheme presented 
in Section III-DI where we sample every T seconds. Upon 
each sample we reset the integrator, so that the mth sample 
correspond to an integral over the interval The resulting 
sampling scheme is depicted in Fig. [3] Since the mth sample is 
influenced by the interval /,„ only, the infinite problem may 
be reduced into a sequence of finite streams of pulses. The 
resulting samples are given by 



c[to] = Sx[m], 



(14) 



where the vector x[to] contains the Fourier series coefficients 
of the signal within the mth interval, Im- As long as S is 
chosen so that it is left invertible, we can obtain the sequence 
of Fourier series coefficients by x[to] = S^c[to]. Extending 
([7|i to the infinite case we obtain: 



y[TO] = H "'"x[to] = V(t[TO])a[m], 



(15) 
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t = mT 



Ci [m] 



x{t)- 



keic 



keic 



T 

]'¥kt 



t = mT 



Fig. 3. Extended sampling scheme using modulating waveforms for an 
infinite pulse stream. 



where t[m] and a[m] are the times and amphtudes of the 
pulses in the interval /„, respectively, and the matrix V 
remains as in (|6]l. For each m, ( fTSl l represents a sum of 
cisoids problem, and thus may be solved as long as K > 2L. 
By choosing p = K = 2L we present a sampling scheme 
which operates at the rate of innovation, and allows for perfect 
reconstruction of an infinite stream of pulses. 
We state our result in a theorem. 

Theorem 2. Consider an infinite stream of pulses given by 
x{t) ^^aih{t~ti), tjGR, azeC. 

where h{t) is a known pulse shape. Assume that there 
are no more than L pulses within any interval I„i — 
[(m — l)r, mT] , m G Z, and that condition (O holds for all 
intervals Im, m £ Z. Choose a set fC of consecutive indices 
for which H{2'Kk/T) 0, Vfc 6 IC. Consider the multichannel 
sampling scheme depicted in Fig. \3\ for some choice of 
coefficients {sik}keK, i = 1, ■ . ■ ,p. Then, the signal x{t) can 
be perfectly reconstructed from the samples {ci[m\}^^-^ with 



T 



x{t) 



E 

fce/c 



'-kt 



dt, ni e Z, (16) 



as long as p > |/C| > 2L, and the coefficients matrix S in 
(II 11 1 is left invertible. 

To the best of our knowledge Theorem |2] presents the 
first sampling scheme for pulse streams with arbitrary shape, 
operating at the rate of innovation. Furthermore, as we show 
in simulations, our method is more stable than previous 
approaches. 



B. Stream of Pulses with Shift-Invariant Structure 

We now focus on a special case of the infinite model ( ITsT l. 
proposed in [21], where the signal has an additional shift- 
invariant (SI) structure. This structure is expressed by the fact 
that in each period T, the delays are constant relative to the 
beginning of the period. Such signals can be described as 



x{t) 



L 



2^ 2_^ai[m]h{t- 



U~mT),te£l CZ[0,T), (17) 



where ai[n] G £2 denotes the £th pulse amplitude on the mth 
period. Assuming condition ^ holds here as well, dTsl l can 
be rewritten as 



y[m] = V(t)a[m], 



(18) 



since now the relative delays in each period are constant. Here, 
a[m] denotes the length-i vector with £th element a^[m]. 

Clearly, the condition for the general model p > 2L is a 
sufficient condition here also, however the additional prior on 
the signal's structure can be used to reduce the number of 
sampling channels. The results obtained in II2TI . for a similar 
set of equations, provide the following sufficient condition for 
unique recovery of the delays and vectors a[m] from ( fTSl i: 



K >2L-T] + 1, 



where 



T] ~ dim (span ({a[; 



GZ})) 



(19) 



(20) 



denotes the dimension of the minimal subspace containing the 
vector set {a[TO], m G Z}. This condition implies that in some 
cases K, and eventually the number of channels p (since p > 
K), can be reduced beyond the lower limit 2L for the general 
model, depending on the value of rj. 

Similar to |21|, recovery of the delays from ( fTSl l can be 
performed using the ESPRIT 128 1 or MUSIC |29| algorithms. 
These approaches, known as subspace methods, require that 
f] — L. In this case they achieve the lower bound of ( fT9] l, 
namely recover the delays using only p > L + 1 sampling 
channels. In cases where r/ < L, an additional smoothing [301 
stage is required prior to using the subspace methods, and 
p > 2L sampling channels are needed. 

To conclude, when the pulse amplitudes vary sufficiently 
from period to period, which is expressed by the condition 
rj — L, the common information about the delays can be uti- 
lized to reduce the sampling rate to {L + 1)/T. Moreover, the 
approach presented here can improve the delays estimation in 
the presence of noise, compared to the one used for the general 
model, since it uses the mutual information between periods, 
rather than recovering the delays for each period separately. 
Will demonstrate this improvement in Section IVI-BI 

This result is summarized in the following theorem. 

Theorem 3. Consider the setup of Theorem |2] where now 

L 

x{t) ai[m]h{t -U- mT), ti e I C [0, T). 

mez e=i 

The signal x(t) can be perfectly reconstructed from the sam- 
ples {ci[m]}^^^, m G Z fl^ long as the coefficients matrix S 
in (II 11 1 is left invertible and 



when rj = L 
when rj < L, 



where ij — dim {span {{a.[m\,m G Z})) denotes the dimen- 
sion of the minimal subspace containing the vector set 
{a[TO], m G Z}. 
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C. Channel Synchronization 

The sampling scheme of Fig. |3] has two main disadvan- 
tages relative to single channel-based schemes: each sampling 
channel requires additional hardware components, and precise 
synchronization of the channel's sampling times is required. In 
this subsection we treat the synchronization issue, and discuss 
the approaches to overcome it. 

One way to synchronize the channels, is on the hardware 
level, for example by using a zero-delay synchronization 
device ll24l . ||31i . Such a device produces accurate trigger 
signals for the samplers and integrators in all the channels. An 
alternative approach is to perform a prior calibration process, 
in which the relative delay of each channel is measured. The 
calibration can be performed at the system manufacturing 
stage or during its power-on, by stimulating the system with a 
known signal. As we now show, once the time offsets between 
the channels are known, they can be compensated. 

Suppose that the ith channel has a time offset of A; £ 
[— Amax, Amax] relative to the optimal sampling instants t ~ 
nT, where A^ax is the maximal possible offset. We assume 
that for each time interval T 

/^(t - to = 0, Vt ^ [A„ax, T - A,„ax) l = l...L. (21) 

This condition ensures that the intervals can be processed 
independently using the approach we now propose. 

The effective delay of the Ith pulse measured in the ith 
channel is ti — A^. Substituting into the fcth Fourier 
coefficient measured in the ith channel satisfies 



T V T 



aie 



j^k{ti-A,) 



1=1 



(22) 



Therefore, from (fTOl i, 

c^ = Y, s^fee-^'^'=^'X[fc] = S^kX[k], (23) 



where we defined 



^ik — ^ik^ 



(24) 



From (|24] | we conclude that when the time offsets between 
the channels are known, the effective mixing matrix of the 
system is S, a matrix whose ifcth element is s^fc. The mis- 
alignment between the sampling channels can be compensated 
by inverting S. We will further discuss the effects of channel 
misahgnment, with unknown delays, in Section IVI-CI 



IV. Modulation Waveforms 



We now thoroughly treat the example in Section III-DI of 
generating the modulation waveforms using periodic signals. 
We first address the case of general periodic waveforms, and 
then focus on the special case of pulse sequences. 



A. General Periodic Waveforms 

Our aim is to show how to obtain the required modulating 
waveforms ^ using a set of p periodic functions, given by 
Pi{t). Such waveforms can be expressed using their Fourier 
series expansion as 



(25) 



fcGZ 



where the fcth Fourier series coefficient of Pi{t) is given by 

dr[k] = - / p,it)e-^^'''dt. (26) 
^ Jo 



The sum in (1251 1 is generally infinite, in contrast to the finite 
sum in (|9]l. Therefore, we propose filtering Pi{t) with a filter 
g{t) which rejects the unwanted elements in the sum ( l25T l. The 
filtered waveforms at the output of g{t) are given by 

Mt)^p^{t)*g{t), (27) 

and are also periodic. Therefore they can be written as 

Mt) = Y,Mky'-^''\ (28) 

feez 



where it can be easily verified that 

l[k] = d,[k]-G(^k 



(29) 



Here G{uj) denotes the CTFT of g{t). From ( l29l l. the shaping 
filter g{t) has to satisfy 

{nonzero uj = ^k, fc € /C 
uj = ^k, k^iC (30) 

arbitrary elsewhere, 

so that di[k] = for k ^ K.. This condition is similar to that 
obtained in |fT3| for single channel sampling. Therefore, the 
class of filters developed there, can also be used here as a 
shaping filter. 

Note that ( [30] l implies that the frequency response of g{t) 
is specified only on the set of discrete points ^k, fc G Z, 
offering large freedom when designing a practical analog filter. 
For instance, when implementing a lowpass filter (LPF) this 
allows a smooth transition band between the passband and the 
stopband of the filter, with a width of 

The resulting scheme is depicted in Fig.|4] The correspond- 



x{t)- 



pi{t)^ g{t) 




■ci 



Pj,{t) g{t) 

Fig. 4. Proposed sampling sclieme, using modulating waveforms. 
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ing elements of the mixing matrix S are given by 

S,k=d,[-{k-[K/2\)]. 



(31) 



The invertibility of S can be ensured by proper selection of the 
periodic waveforms Pi{t). In the next subsection we discuss 
one special case, which allows simple design S. 

B. Pulse Sequence Modulation 

We follow practical modulation implementation ideas pre- 



sented in I22i . |24J, and consider the set of waveforms 

N-l 
mGZ n— 

for i — 1, . . . where p{t) is some pulse shape and ai[n\ 
is a length-A^ sequence. Our aim is to calculate the mixing 
matrix S, when using the filtered version of ( |32] | as modulating 
waveforms. To this end, we first compute the Fourier series 
coefficients di[k] of Pi{t) as 

N-l T 
N-l ^-(m-l)T 

„=o mez"'-™'^ 

^ a,[n] / p{t - nT/N)e-^'^^'dt 

n=0 

=l£'a,,Mpf^fc)e-^^'^", (33) 



n=0 



where denotes the CTFT of p{t). 

Combining (|33] | with (1311 1 and 



n=0 \ / \ / 

where we defined k' = — (fc — [K/2\). The resulting matrix 
S can be decomposed as 



S = AW*, 



(35) 



where A is a p x matrix with mth element equal to ajn], 
W is an X iiT matrix with nkth element equal to e~^^^ ", 
and $ is a X if diagonal matrix with fcth diagonal element 



(36) 



From this decomposition it is clear that $ has to be invertible 
and the matrix AW has to be left invertible, in order to 
guarantee the left invertibility of S. We now examine each 
one of these matrices. 

We start with the matrix From ((301), G {^k) ^ for 
k € IC. Therefore, we only need to require that P {^k^ ^ 
for fc G in order for $ to be invertible. A necessary 
condition for the matrix AW to be left invertible, is that W 
has full column rank. The matrix W is a Vandermonde matrix, 
and therefore has full column rank as long as iV > if ll32ll . 
The left invertibility of the px K matrix AW can be ensured. 



by proper selection of the sequences a^JT^], where a necessary 
condition is that p > K. 

We summarize our results in the following proposition: 

Proposition 1. Consider the system depicted in Fig. |4] where 
the modulation waveforms are given by ( I32I ). If the following 
conditions hold 

1) p>K,N> K, 

2) The frequency response of the shaping pulse g{t) satis- 
fies 

3) The frequency response of the pulse p{t) satisfies 
P{^k) ^ Ofor keIC, 

4) The sequences ai[n\ are chosen such that the matrix 
AW has a full column rank, 

then the mixing matrix S in ( 1111 ) is left invertible. 

We now give two useful configurations, that satisfy the 
conditions of Proposition [T] 

7 ) Single Generator: We create the ith sequence, by taking 
a cyclic shift of one common sequence a[n\ as 



[n] = a[n — i + 1 mod N], 



(37) 



where we assume p — N. Clearly, the corresponding wave- 
forms can be created by using only one pulse generator, 
where the waveform at the ith channel is a delayed version 
of the generator output, delayed by {i — 1)T/N time units. 
This suggests, that in contrast to the direct scheme in Fig. [T] 
which requires multiple frequency sources, here only one pulse 
generator is required which simplifies the hardware design. It 
is easy to see that with this choice, A will be a circulant 
matrix. Such a matrix can be decomposed ll33l as 



A = F^diag (Fa) F, 



(38) 



where F is a iV x TV unitary discrete Fourier transform (DFT) 
matrix, and a is a length-iV vector containing the elements of 
the sequence a[n]. Therefore, for A to be invertible the DFT 
of the sequence a[n] can not take on the value zero. 

We now give an example for such a selection of the system's 
parameters. We set p ~ N ~ K, and choose 



The frequency response of this pulse satisfies 

T 

p J 2N'^ . sine I - 

N 



T 

Piuj) = —e 

^ ' T\T 



2'kN 



Therefore, 





T . 


f k\ 




= — smc 






N 





(39) 



(40) 



(41) 



which is non-zero for /c G /C. In addition we choose the 
sequences ai\v\ as sequences of ±ls, created from cyclic 
shifts of one basic sequence, in a way that yields an invertible 
matrix A. Such rectangular pulses with alternating signs can 
be easily implemented in hardware [22 1. In Figs. |5] and |6] 
one modulating waveform is shown in the time and frequency 
domains, for p — N ^ K ~ 7. The original time-domain 
waveform is comprised of rectangular pulses, whereas lowpass 



g 



filtering results in a smooth modulating waveform. Switching 
to the frequency domain, the Fourier series coefficients are 
shaped by P{uj), the CTFT of the pulse shape. The shaping 
filter frequency response, G{uj), is designed to transfer only 
the Fourier coefficients whose index is a member of the set 
/C = {—3, . . . , 3}, suppressing all other coefficients. 
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Fig. 5. Modulating waveform in tlie time domain, before and after filtering. 



from the previous discussion, we demonstrate how to reduce 
the number of required generators, for the current setting. 
For simplicity, we assume p > 2pe and that two different 
generators are used. The first half of the sampling channels use 
delayed versions of the first generator output, and the second 
half uses the second generator By proper selection of the two 
sequences, the condition mentioned above can be satisfied. 

We now give a numerical example for a such choice. 
We assume i = 4 and set = /v = 9 and p = 18 
sampling channels, which are based on two generators only. 
Each generator produces a different sequence of ±1, chosen 
randomly. In Fig. |7] we plot the log of the maximal condition 
number of S, obtained when going over all possible options 
for omitting pe rows from S. It can be seen that for Pe < 6 
a relatively low condition number of the matrix S is achieved 
in the worst case, suggesting that its rank is N as required 
and is not ill-conditioned. Therefore, we can overcome 
failure in up to 6 sampling channels, using this system. In 
this case the required K ~% consecutive Fourier coefficients 
can be obtained from the remaining channels, allowing the 
perfect recovery of the L — A pulses. In contrast, when 
channels fail in the direct scheme of Fig. [T] a set of ii' = 9 
consecutive Fourier coefficients cannot always be obtained. 
Therefore, perfect recovery of the signal using methods such 
as annihilating filter or matrix pencil is not guaranteed. 




-5 5 

frequency [2%/J] 

Fig. 6. Modulating waveform in the frequency domain. 
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Fig. 7. Robustness to sampling channels failure example, N = K 
p = 18, two pulse generators. 



2) Robustness to Sampling Channels Failure: Next we 
provide a setup which can overcome failures in a given 
number of the sampling channels. The identification of the 
malfunctioning channels is assumed to be performed by some 
external hardware. 

We consider a maximal number of malfunctioning chan- 
nels and assume that p > N + Pf.- Ivi order to ensure unique 
recovery of x from the remaining channels, the submatrix 
S, obtained from S by omitting of the corresponding rows, 
should be left invertible. Since this has to be satisfied for every 
possible selection of p^ rows, we need to design S such that 
any p—pe rows will form a rank-X matrix. Following our ideas 



V. Related Work 

A. Single-Channel Sampling with the SoS Filter 

The work in 1131 considered single-channel sampling 
schemes for pulse streams, based on a filter which is comprised 
of a Sum of Sines (SoS) in the frequency domain. This filter 
can be expressed in the time domain as 

5(t) = rect('|W6fee^^'=S (42) 

where /C is the chosen index set, and the coefficients {hk]k(^K 
have arbitrary nonzero values. 
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To explore the relation of our method to |fT3]| . we first focus 
on periodic streams of pulses with period T. By sampling a 
periodic stream of pulses using the scheme depicted in Fig. [8] 
the following samples are obtained 



c n 



fce/c 



fcriTa 



0, 



(43) 



where Tg — T/p is the sampling period. Using the matrix V 
defined in (|6]l only now with parameter ts = {0, Ts, . . . , (p — 
l)Ts}, and defining the diagonal matrix B with kt\\ diagonal 
element bk, ( |43] | can be written in matrix form as 



c = V(-ts)Bx. 



(44) 



Therefore, this is a special case of our multichannel sampling 
scheme in ( fTTI ) with mixing matrix S = V(— ts)B. The matrix 
B is invertible by construction, and V(— t^) is a Vandermonde 
matrix with distinct times, so that S is left-invertible as long 
as p > if. Using this choice of S, the p samples taken over 
one period in ||T3]| . are equal to the samples at the output of 
the p channels in our scheme. 



x{t) 



SoS Filter 



■ c n 



t = nT, 



Fig. 8. Single-channel sampling scheme using the SoS filter. 

Exploiting the compact support of the SoS filter, the method 
was extended to the finite and infinite settings as well 1 13|. The 
extension is based on using an r-fold periodic continuation 
of the SoS filter g{t), where the parameter r depends on 
the support of the pulse-shape h{t). For the infinite case, the 
samples of the two schemes still coincide. However, since ITSl 
is a filtering based scheme, proper separation of at least 1.5T 
between periods is required in order to obtain independent 
processing of each period. While keeping the same sampling 
rate of 2L/T, the model in lfT3l has a rate of innovation of 
2L/(2.5T). Hence the single channel configuration of lT3l 
does not achieve the rate of innovation in the infinite case. 

We now examine the modulation waveforms that result from 
ifTSJ . Since S = V(-t,)B, 



fce/c 



-3^k{t-iT/p) 



(45) 



It is easily shown that these waveforms can be expressed as 



s,it)^~g{~it~iT/p)), 



(46) 



where g{t) is the periodic continuation of the SoS filter g{t). 
Therefore, in each channel the signal is modulated by a 
delayed version of the periodic SoS filter The equivalence 
of the schemes is easy to explain: sampling the convolution 
between the input signal and the SoS filter in 1 13|, is equivalent 
to performing inner products (multiplication followed by inte- 
gration) with delayed and reflected versions of this filter This 
relation provides another valid class of modulation waveforms. 



B. Multichannel Schemes for Shift-Invariant Pulse Streams 

Another related work is |21| which treats the SI signal 
model ([TtT i presented in Section IIII-BI The sampling scheme 
proposed in ||2T1 is depicted in Fig. |9] In each channel, the 
input signal is filtered by a band-limited sampling kernel 
s|(— t) followed by a uniform sampler operating at a rate 
of l/T. After sampling, a properly designed digital filter 
correction bank, whose frequency response in the DTFT do- 
main is denoted here by 'M.{e''^'^), is applied on the sampling 
sequences. The exact form of this filter bank is detailed in 
[21 J. It was shown in [21], that the ESPRIT algorithm can 
be applied on the corrected samples, in order to recover the 
unknown delays. 



4 i-t) 



( = riT 



t = nT 



■ch[n] 
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Fig. 9. Proposed sampling scheme in 1211 . 

The sampling rate achieved by the method in 1*211 is 
generally 2L/T, where for certain signals it can be reduced to 
(i + l)/T. Such signals satisfy dim(span({d[m]. m e Z})) = 
L, where the vectors d[m] contain the samples at the output of 
the scheme depicted in Fig. |9l This condition is different than 
iT% . which directly depends on the vectors a[m]. Therefore 
the sampling rate, when using the scheme in (2T\, can be 
reduced to {L + 1)/T for different signals. This fact is not 
surprising, since each approach has a different analog sampling 
stage. In both methods, the worst-case minimal sampling rate 
is 2L/T. 

The approach in 11211 has two main advantages over the 
proposed method in this work. The first is that condition (|2]l is 
not required. Therefore ||2T| can also treat pulses with infinite 
time support, in contrast to our method. Another advantage is 
that it can support single channel configurations. It was shown 
in 11271 . ll34l . that one sampling channel followed by a serial to 
parallel converter, can be used in order to produce the parallel 
sampling sequences in Fig.|9] 

On the other hand, the method depicted in Fig. [3] has several 
advantages over 11211 . First, the equivalent stage for the digital 
correction in fTP\, is replaced by inversion of the matrices H 
and S, since 



y[m] = (SH)-icH- 



(47) 



This operation can be viewed as a one-tap digital correction 
filter bank, in contrast to the filter M(eJ"-'^), which generally 
has a larger number of taps. Therefore, the proposed correction 
stage, is much simpler and requires lower computational 
complexity, than the one in [21 j. 

An additional advantage of our scheme, is that the approach 
of f2V\ requires collection of an infinite number of samples, 
even when the input signal contains a finite number of periods. 
This requirement is due to the infinite time support of the band- 
limited sampling kernels. Moreover, if one is interested only 
in a finite time interval of the signal, the method in 1211 does 
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not allow processing it separately. This is in contrast to the 
proposed scheme, which integrates finite time intervals, and 
can collect samples only from the relevant periods. We will 
demonstrate this advantage in Section IVII 



C. Modulated Wideband Converter 

The concept of using modulation waveforms, is based on 
ideas which were presented in f22l-f24l, f35l. We now briefly 
review the sampling problem treated in [241 and its relation to 
our setup. As we show the practical hardware implementation 
of both systems is similar. 

The model in [24] is of multiband signals: signals whose 
CTFT is concentrated on A'bands frequency bands, and the 
width of each band is no greater than B. The location of the 
bands is unknown in advance. An example of such a signal 
is depicted in Fig. [10] A low rate sampling scheme allowing 
recovery of such signals at a rate of 4_BA^bands was proposed in 
Q. This scheme exploits the sparsity of multiband signals in 
the frequency domain, to reduce the sampling rate well below 
the Nyquist rate. In Il22l - ll24l . this approach was extended to 




Fig. 10. Multiband signal model. 

a more practical sampling scheme, which uses a modulation 
stage and referred to as the Modulated Wideband Converter 
(MWC). In each channel of the MWC, the input is modulated 
with some periodic waveform, and then sampled using a LPF 
followed by a low rate uniform sampler. The main idea is that 
in each channel, the spectrum of the signal is scrambled, such 
that a portion of the energy of all bands appears at baseband. 
Mixing of the frequency bands in 1241 is analogous to mixing 
the Fourier coefficients in Fig. [3] 

We note here some differences between the methods. First, 
following the mixing stage, we use an integrator in contrast 
to the LPF used in [24J. This difference is a result of the 
different signal quantities measured: Fourier coefficients in our 
work as opposed to the frequency bands content in |24|. The 
second difference is in the purpose of the mixing procedure. 
In Il24l mixing is performed in order to reduce the sampling 
rate relative to the Nyquist rate. In our setting, the mixing is 
used in order to simplify the hardware implementation and to 
improve robustness to failure in one of the sampling channels. 

Nonetheless, the hardware considerations in the mixing 
stage in both systems is similar Recently, a prototype of the 
MWC has been implemented in hardware [22|. This design is 
composed of p = 4 sampling channels, where the repetition 
rate of the modulating waveforms is 1 /T « 20 MHz. In each 
period there are N = 108 rectangular pulses. This prototype, 
with certain modifications, can be used to implement our 
sampling scheme as well. These modifications mainly include 
adding shaping filters on modulating waveforms lines, and 
reducing the number of rectangular pulses in each period. 



VI. Simulations 

In this section we provide several experiments in which we 
examine various aspects of our method. The simulations are 
divided into 4 parts: 

1) Evaluation of the performance in the presence of noise, 
and comparison to other techniques; 

2) Demonstration of the recovery method for pulse streams 
with SI structure; 

3) Evaluation of the effects of synchronization errors be- 
tween the channels; 

4) Examination of the use of practical shaping filters. 

A. Performance in the Presence of Noise 

We demonstrate the performance of our approach in the 
presence of white gaussian noise, when working at the rate of 
innovation. We compare our results to those achieved by the 
integrators |19| and exponential filters |20| based methods, 
since these are the only approaches which can work at the 
same rate, for infinite stream of pulses. 

We examine three modulation waveforms presented in Sec- 
tions UPD] and EA] cosine and sine waveform (tones), filtered 
rectangular alternating pulses (rectangular) and waveforms 
obtained from delayed versions of the SoS filter (SoS). For the 
rectangular pulses scheme, the modulation waveforms are gen- 
erated using a single generator, as discussed in Section HV-BI 
The shaping filter g{t) is an ideal LPF with transition band 
of width 2-k/T. Following [l20l, the parameters defining the 
impulse response of the exponential filters are chosen as 
a = 0.2T and /? = 0.8T. 

We focus on one period of the input signal, which consists 
of L = 2 Diracs with t = [0.256r, O.SST]'^, and amplitudes 
a — [1, 0.8]"^. We set the system parameters as p ~ K = N = 
5. The estimation error of the time-delays versus the SNR is 
depicted in Fig. [TT] for the various approaches. Evidently, 
our technique outperforms the integrators and exponential 
filters based methods in terms of noise robustness, for all 
configurations. There is a slight advantage of 2dB for the 
schemes based on tones and SoS, over alternating pulses, 
where the first two configurations have similar performance. 

Turning to higher order problems, in Fig. [12] we show the 
results for L = 10 Diracs with times chosen in the interval 
[0,r) and amplitudes equal one, with N = p ^ K — 21. 
The instability of the integrators and exponential filters based 
methods becomes apparent in this simulation. Our approach in 
contrast achieves good estimation results, demonstrating that 
our method is stable even for high model orders. 

The performance advantage of the tones and SoS based 
schemes is now around 3.5dB. We conclude that from a noise 
robustness point of view, using multiple frequency sources or 
SoS waveforms is preferable over a single pulse generator. 
However, as discussed in Section III-DI pulse sequences based 
schemes can be advantageous from practical implementation 
considerations, and reduce the hardware complexity. In addi- 
tion, the performance degradation is reasonable, and the esti- 
mation error is still significantly lower than that of competing 
approaches. Therefore, the flexibility of our architecture, al- 
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Fig. 1 1 . Performance in tlie presence of noise, at the rate of innovation. Tlie 
signal consists of L = 2 Diracs. 



lows the system designer to decide between better performance 
in the presence of noise, or lower hardware complexity. 
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Fig. 12. Performance in tlie presence of noise, at tlie rate of innovation. The 
signal consists of L = 10 Diracs. 

Next, we compare our scheme to the one presented in lfT6ll . 
This approach, which is based on B-splines lf36l and E-splines 
(J? I sampling kernels, operates at a rate higher than the rate of 
innovation. According to the main theorem in lfT6l . an infinite 
stream of Diracs is uniquely determined from uniform samples 
taken at the output of a B-spline or E-spline sampling kernel, 
if there are at most L Diracs in an interval of size 2LSts- Here 
is is the sampling interval, and S" = 2L is the time support 
of the sampling kernel. In our setting there are L Diracs in 
an interval of size T requiring at least p = (2L)^ samples per 
period T . 

We choose L = 4 Diracs, with delays t = 
[0.213T, 0.452r, 0.664T, 0.745r]^ and amplitudes a = 
[1,0.9,0.7,0.6]"^. We compare our tones based configuration 



to the both B-spline and E-spline techniques of lfT6l . The 
parameters defining the E-spline kernel fT6l, fJTl were chosen 
in order to obtain real valued sampling kernels, and were tuned 
empirically to obtain the best performance. For all algorithms 
p = 64 samples are used. In order to exploit the oversampling 
in our approach, the Kumaresan and Tufts method 1,27 J is used 
for the delays recovery. The estimation error of the time-delays 
versus SNR is depicted in Fig. [13] From the figure it can be 
seen that our scheme exhibits better noise robustness than both 
B-spline and E-spline based methods. 




40 60 
SNR [dB] 

Fig. 13. Performance in the presence of noise, using p = 64 samples. The 
signal consists of L = 4 Diracs. 



B. Sampling of Pulses with SI Structure 

We now consider sampling of streams of Diracs with SI 
structure. We compare our method to the one presented in 
lim . For our scheme, we examine two recovery options. The 
first is to process each period separately, namely, to recover 
the delays from each period independently (standard recovery). 
The second is to follow the approach presented in Section lTlI-BI 
and to recover the common delays from all periods using the 
ESPRIT [28 1 algorithm (SI recovery). 

We consider 25 periods with L = 4 Diracs per period, 
and relative delays of t = [0.213T, 0.452T, 0.664r, 0.745r]^. 
The amplitudes in each period were taken as an inde- 
pendent Gaussian random variables, with means fi^ = 
[1,0.9,0.7,0.6]^ and standard deviation a = 0.1 For the 
method in ||2T1 . we chose a single channel scheme with an 
ideal LPF as sampling kernel. For both methods p = 8 
samples per period were taken. The technique in [21) requires 
theoretically infinite number of samples, due to the infinite 
time support of the sampling kernel. However, to compare 
between the two methods, we used only 8 ■ 25 samples (25 
periods, with 8 samples per period). 

The estimation error of the time-delays versus the SNR 
is depicted in Fig. [141 For SNR levels above 15dB, there 
is a clear advantage to the SI recovery method, over the 
standard approach. Hence, as expected, the use of the mutual 
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information between periods on the delays, improves the 
estimation significantly. Comparing the performance of our 
approach to the one of |21|, it can be seen that up to SNR 
levels of 20dB both methods achieve similar performance. 
However, for higher SNRs, the method in fIT] suffers from 
a dominant error caused by the fact that only a finite number 
of samples were used. This demonstrates the advantage of our 
scheme, which operates on finite time intervals, in cases where 
the signal consists of a finite number of periods. 
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Fig. 14. Sampling of pulses with SI structure. The signal consists of 25 
periods, each with L = 4 Diracs. 



C. Synchronization Errors 

We now study the error caused by synchronization errors 
between the channels. We consider sampling of L = 4 
Diracs, with delays t = [0.213T, 0.452r, 0.664T, 0.745T]^ 
and amplitudes a ~ [10.90.70.6]"^, using p = 9 sampling 
channels. We set the sampling time of the last channel to be 
shifted by An,ax relative to the first channel. The offsets of the 
other channels are drawn uniformly on the interval [0, A^ax]- 
We plot the standard deviation of the estimated time delays 
error, as a function of the maximal offset Amax, for different 
SNR values. We note that in this experiment, the misalignment 
between the sampling channels is not being compensated using 
the approach discussed in Section IIII-CI 

The results are shown in Fig. [15] The dash-dotted line 
denotes the linear curve y = x. First, it can be seen that 
when the synchronization error is less than 10 percent of the 
estimation error, the synchronization error is negligible, and 
the error is mainly due to the noise. When the synchronization 
error becomes large, the time delay estimation error degrades 
linearly. In general, the estimation error is bounded from below 
by the synchronization error (for A^ax < 0.03T). 
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Fig. 15. Standard deviation of the time-delays estimation error, as a function 
of the maximal offset between channels, for various SNR levels. 



coefficients whose index is not in the set JC is not perfect. 
We set the shaping filter g{t) to be a Chebyshev (Type I) 
LPF ll38l of various orders, with ripple 3 dB. The Chebyshev 
filter is a good choice for our requirements since it has a 
steeper roll-off than other filters, resulting in better rejection 
of the undesired coefficients. The rapid transition between the 
pass-band and stop-band of the Chebyshev filter comes at the 
expense of larger ripple in the pass-band, however, ripple is 
of minor concern for our method since it is digitally corrected 
when inverting the matrix S. The cutoff frequency was set 
to ^lK/2\. The frequency response of the various filters is 
shown in Fig. [16] 



d. 

ideal filter 
Chebyshev order 3 
Chebyshev order 10 




frequency [27t/T] 

Fig. 16. Frequency responses of the shaping filters: ideal shaping filter vs. 
practical Chebyshev filters. 



D. Practical Shaping Filters 

We now explore the use of practical shaping filters, for 
the rectangular pulses scheme, rather than the ideal ones 
used above. Once practical filters are used, the rejection of 



The estimation error of the time-delays versus the SNR is 
depicted in Fig. [17] for various filter orders. The simulation 
consists of L = 2 Diracs with t = [0.256T, 0.46T]^, and 
amplitudes a = [1,0.8]^. Clearly, a Chebyshev filter of order 
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10 closely approaches the performance of an ideal LPF. In 
addition, for SNR levels below 50dB, using a Chebyshev 
filter of order 6 provides good approximation. Therefore, 
the modulation waveform generation stage of our proposed 
method can be implemented using practical analog filters. 
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Fig. 17. Performance of practical shaping filters of various orders vs. ideal 
filtering. 



VII. Conclusion 

In this work, we proposed a new class of sampling schemes 
for pulse streams. Our approach allows recovery of the delays 
and amplitudes defining such a signal, while operating at the 
rate of innovation. In contrast to previous works lfT8l - ll20l 
which achieved the rate of innovation, our approach supports 
general pulse shapes, rather than Diracs only. In addition, as 
we demonstrate by simulations, our method exhibits better 
noise robustness than previous methods lfT6l . lfT9l . ||20I . and 
can accommodate high rates of innovation. 

The proposed scheme is based on multiple channels, each 
comprised of mixing with a properly chosen waveform fol- 
lowed by an integrator We exploit the degrees of freedom in 
the waveforms selection, and provide several useful configu- 
rations, which allow simplified hardware implementation and 
robustness to channel failure. Using simulations we further 
explored practical issues, such as effects of misalignment 
between the sampling channels and usage of standard analog 
filters, in the waveform generation stage. 

Our method can be viewed as a part of a broader framework 
for sub-Nyquist sampling of analog signals, referred to as 
Xampling 1221 . ||23]| . We draw connections with the work 
in Il22ll - ll24l . which proposed a Xampling architecture for 
multiband signals. We showed that the hardware prototype of 
the analog front-end, implemented for the multiband model, 
can be used in our scheme as well with certain modifications. 
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